Persistence in the Zero-Temperature Dynamics of the Random Ising Ferromagnet on 

a Voronoi-Delaunay lattice 
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The zero-temperature Glauber dynamic is used to investigate the persistence probability P(t) 
in the randomic two-dimensional ferromagnetic Ising model on a Voronoi-Delaunay tessellation. 
We consider the coupling factor J varying with the distance r between the first neighbors to be 
J(r) oc e~ ar , with a > 0. The persistence probability of spins flip, that does not depends on time 
t, is found to decay to a non-zero value P(oo) depending on the parameter a. Nevertheless, the 
quantity p(t) = P(t) — P(oo) decays exponentially to zero over long times. Furthermore, the fraction 
of spins that do not change at a time t is a monotonically increasing function of the parameter a. 
Our results are consistent with the ones obtained for the diluted ferromagnetic Ising model on a 
square lattice. 

PACS number(s): 05.20-y, 05.50+q, 64.60.Cn, 75.10.Hk, 75.40.Mg 
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I. INTRODUCTION 

The Lenz-Ising model is probably the simplest non- 
trivial model for cooperative behavior that shows sponta- 
neous breaking of symmetry . Because of this simplicity 
and the fact that each individual element of the model 
modifies its behavior according to the other individuals in 
its vicinity, it has a vast number of applications ranging 
from solid-state physics to biology [1-3]. An important 
extension to this model is the introduction of disorder 
such as; random external fields, random exchange pa- 
rameters, and dilution, where only a fraction p of the 
lattice sites are occupied by spins. In particular, the two 
dimensional diluted ferromagnetic random Ising model 
is very important to describe magnetic properties of sev- 
eral condensed matter systems. Moreover, it constitutes 
a marginal situation of the Harris criterion [4], and for 
many years there has been several theoretical and numer- 
ical studies [5-8] in order to clarify the properties of the 
dilution model. 

On the other hand, only in the last ten years the "per- 
sistence" problem in the randomic two-dimensional ferro- 
magnetic Ising model has attracted considerable interest 
[9-11]. In this most general form, this problem involves 
the fraction of space which persists in its initial state 
until some time later. Hence, in the non-equilibrium dy- 
namics of spin systems we are interested in the fraction 
of spins P(t), that persist in the same state as at t = 
up to some later time t. For the pure ferromagnetic d- 
dimensional (d < 4) Ising model, P(t) has been found to 
decay algebraically [9,12] 



P(t) cx t 



-9{d) 



(1) 



where 9(d) is a non-trivial persistence exponent. For the 
same model [10] at higher dimensions (d > 4) and for the 
two-dimensional ferromagnetic g-state Potts model [13] 
(q > 4) it has been shown through computer simulations 



that P(t) is not equal to zero for t — > oo. This charac- 
terictic is some time called as "blocking". Therefore, if 
P(oo) > we can reformulate the problem by restrict- 
ing our observation to those spins that eventually flip. 
Hence, we can consider the behavior of 



p(t) = P(t) - P(oo) 



(2) 



By considering the dynamics of the local order param- 
eter the persistence problem can be generalized to non- 
zero temperatures [14,15] Recently [16,17] the attention 
has been turned to the persistence problem in systems 
containing disorder. Numerical simulations of the zero- 
temperature dynamics of the bond diluted (weak dilu- 
tion [16] or strong dilution [18]) two-dimensional Ising 
model also reported "blocking" evidences. Howard [19] 
has found evidence of an exponential decay of the persis- 
tence with blocking for the homogeneous ferromagnetic 
Ising model on the homogeneous tree of degree three 
(T) with random spin configuration at time 0. Here we 
present results of an extensive numerical study of the per- 
sistence of the ferromagnetic Ising model on the Voronoi- 
Delaunay lattice. In this lattice, the coordination number 
and the distance between the first neighbors sites is ran- 
domic. As the bond length between the first neighbors 
varies randomly from neighbor to neighbor, we consider 
that the coupling factor depends on the relative distance 
Tij between sites i and j according to the following ex- 
pression: 



Jij = J e~ 



(3) 



where Jo is a constant and < a < 1. The question to be 
answered here is: does the persistence behavior change 
with this type of randomness or is its behavior the same 
as in the diluted ferromagnetic Ising model in a regular 
lattice? In the present work, we show that, with this 
type of randomness, pit) presents an exponential decay 
as in the strongly diluted ferromagnetic Ising model in 
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two-dimension [17] in contrast with the behavior of the 
pure and weakly diluted models. 



II. MODEL AND SIMULATION 

The Voronoi construction or tessellation for a given 
set of points in the plane is defined as follows. For each 
point, we first determine the polygonal cell consisting of 
the region of space nearer to that point than any other 
point. Whenever two such cells share an edge, they are 
considered to be neighbors. From the Voronoi tessel- 
lation, we can obtain the dual lattice by the following 
procedure: when two cells are neighbors, a link is placed 
between the two points located in the cells. From the 
links, one obtains the triangulation of space that is called 
the Delaunay lattice. The Delaunay lattice is dual to the 
Voronoy tessellation in the sense that points correspond 
to cells, links to edges and triangles to the vertices of the 
Voronoi tessellation. 

We consider a two-dimensional Ising ferromagnetic 
on this Poissonian random lattice which Hamiltonian is 
given by: 



-KH = — JijSiSj 

<i,j> 



(4) 



where Si — ±1 are the Ising spins situated on every site 
of a Delaunay lattice with (LxL — N) sites and periodic 
boundary conditions; K = 1/ksT , T is the tempera- 
ture and Ub is the Boltzmann constant. The summation 
in Eq. (4) runs over all nearest-neighbors pairs of sites 
(points in the Delaunay construction) . In this lattice the 
coordination number varies locally between 3 and oo and 
the coupling factor Jij depends on the distance between 
first neighbors according to Eq. 3. 
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For simplicity, the length of the system is defined here 
in terms of the size of a regular lattice, L — N 1 / 2 . We 
perform simulations over a lattice with L = 500. A 
randomly initial configuration of spins is obtained and 
P{t) is calculated over 6000 of Monte Carlo steps and a 
quenched average is done over 10 different Delaunay lat- 
tices for each Monte Carlo step. The temperature zero 
Glauber dynamics was utilized in order to check the num- 
ber of spins that never change their state at a time t. In 
this dynamics we start with a randomly initial spin con- 
figuration and allow it to be updated by selecting one 
spin to be flipped at random or following a given log- 
ical sequence. The selected spin, Si, is flipped or not 
according to AEi, where AEi, is the energy of site i. If 
AEi < the spin 5, is flipped with probability one. If 
AEi = the spin Si is flipped with probability 1/2 cho- 
sen at random. Finally if AEi > the spin is not flipped. 
One Monte Carlo step corresponds to application of the 
above rule for all spins of the lattice. The system config- 
uration is left to evolve until a given Monte Carlo step 
t = 6, 000. The number, n(t), of sites that do not change 
at this time t is computed for each Monte Carlo step for 
the determination of the persistence probability given by 
[9]: 
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FIG. 2. Fraction of spins that never flip versus a. 
< n(t) > 



P(t) = 



N 



(5) 



FIG. 1. ln[P(t)] versus ln(t) for three different values of 
the parameter a and lattice size L = 500. 



The data have been obtained for 10 a values in the range 
< a < 1. 



III. RESULTS AND CONCLUSION 

For a varying from (where the coupling constant is 
independent of the distance between first neighbors) to 
1 the persistence P(t) seems not to decay algebraically 
before the "freezing" as we can see in Fig. 1. In this 
figure we plot In P(t) versus t for a = , 0.5 and 1. We 
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can also verify that, P(t) = P(oo) for t > t*, where t*(a) 
depends on the parameter a, growing with the a value. 
The fraction of frozen sites (i.e. the sites that never flip) 
in function of the parameter a is shown in Fig. (2). This 
fraction has a maximum value near a = 0.9) decreasing 
for (a = 1). Finally in Fig. (3) we plot In [p(t)] versus 
In t for the same values of a. From this figure we can 
verify that p(t) decay exponentially for long times. This 
result agrees with the results obtained by Newman and 
Stein [17] for the persistence in the strongly diluted Ising 
ferromagnet. This behavior occurs for a = 0. This result 
is contradictory for us once we know that for a = the 
ferromagnetic Ising model on a Dclaunay lattice has the 
same critical exponents that the pure ferromagnetic Ising 
model in a square lattice has [8] . This fact does not agree 
with the persistence behavior of the pure and weakly di- 
luted ferromagnetic Ising model reported by Jain [16]. 
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FIG. 3. A log-log plot of p(t) versus t for three different a. 

In summary we have presented new data for the zero- 
temperature dynamics of the two-dimensional ferromag- 
netic Ising model on a random Poissonian lattice. This 
system exhibit "blocking" what means that the persis- 
tence does not go to zero when t — > oo. We also find that 
p(t) decay exponentially in the long time regime. The 
fraction of "frozen" sites increases from a non-zero value 



with the parameter a. Our results strongly suggest that 
the persistence behavior is not algebraic supporting the 
suggestion that the decay of the persistence probability 
can be not-algebraic for certain classes of models. 
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